Quantum Disordered Ground States in Frustrated Antiferromagnets with Multiple 
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We present a certain class of two-dimensional frustrated quantum Ifeisenberg spin systems with 
multiple ring exchange interactions which are rigorously demonstrated to have quantum disordered 
ground states without magnetic long-range order. The systems considered in this paper are s=l/2 
antiferromagnets on a honeycomb and square lattices, and an s=l antiferromagnet on a triangular 
lattice. We find that for a particular set of parameter values, the ground state is a short-range 
resonating valence bond state or a valence bond crystal state. It is shown that these systems are 
closely related to the quantum dimer model introduced by Rokhsar and Kivelson as an effective 
low-energy theory for valence bond states. 



I. INTRODUCTION 

Quantum frustrated Heisenberg antiferromagnets have 
attracted a great deal of interest in connection with the 
search for exotic phases, such as spin liquids. Ander- 
son proposed a resonating valence bond (RVB) state as 
a prototype of spin liquids three decades agoi. Since 
then, a number of works have studied the conjecture that 
both geometrical frustration and strong quantum fluctua- 
tion may destroy magnetic long-range order and stabilize 

a quantum disordered liquid-like ground state without 
symmetry breaking^i3i4i5,6/,8,9,io,ii,i2,i3,i4,i5,i6,i7_ p^^ 

example, it is believed quite likely that quantum spin 
systems on Kagome lattices and pyrochlore lattices can 
exhibit quantum spin liquid states, though the elu- 
cidation of their ground states remains an important 
unsolved problem. In related works, the quantum 
dimer model (QDM) introduced by Rokhsar and Kivel- 
son has been extensively studied as an effective low- 
energy theory of quantum disordered states with short- 
range antiferromagnetic correlation and a spin excitation 
gg^pi8ji9j2i,22^^ rpj^g important feature of the QDM is 
that its Hilbert space is spanned by only dimer covering 
states composed of singlet pairs of nearest-neighbor spins, 
and this is the key to its success in describing quantum 
disordered states^'*. It was shown in the pioneer paper 
by Rokhsar and Kivelson that a short-range RVB state 
is realized in the QDM as its exact ground state. It was 
also argued by several authors that the QDM may be an 
effective low-energy theory of quantum antiferromagnets 
on the kagome lattice. 

However, it is not unclear how to derive the QDM di- 
rectly from the original quantum spin Hamiltonian by 
truncating the Hilbert space. Indeed, it appears difficult 
to find a general answer to this question. Therefore, in 
this paper, we do not seek such a general result. Instead, 
to obtain better insight regarding the mechanism stabi- 
lizing quantum spin liquids, we attempt to make progress 
toward determining what kinds of quantum spin systems 
possess a low-energy sector described by the QDM. 

The main result of this paper is the rigorous proof that 
a certain class of quantum antiferromagnetic spin systems 



in two dimensions with multiple ring exchange interac- 
tions is equivalent to the QDM. These systems exhibit a 
short-range RVB state or a valence bond crystal (VBC) 
state, depending on the values of their parameters. Our 
results clarify the important role played by the ring ex- 
change interactions in the realization of quantum spin 
liquid states-- . 

This paper is organized as follows. In Sec. II, we present 
results of an s = 1/2 Heisenberg spin system with mul- 
tiple ring exchange interactions defined on a honeycomb 
lattice. It is shown that for a particular parameter the 
rigorous ground state is the short range RVB. In this 
RVB state, spin-spin correlations decay exponentially, 
and therefore there is no magnetic long range order. On 
the other hand, correlations between spin-singlet dimers 
exhibit power law long-distance behavior, implying the 
existence of low-lying spin-singlet gapless excitations. In 
Sec. HI and Sec. IV, we consider an s = 1/2 model on a 
square lattice and an s = 1 model on a triangular lattice, 
respecticely. Discussion and summary are given in Sec.V. 



II. s = 1/2 HEISENBERG ANTIFERROMAGNET 
ON A HONEYCOMB LATTICE 



Model Hamiltonian 



We consider the s = 1/2 c^uantum Heisenberg anti- 
ferromagnet on a honeycomb lattice (see FIG. 1(a)), of 
which the ground state is exactly obtained. The Hamil- 
tonian is given by, H = ffx + -ffR, with 
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Here, (ij), < ij > and <C ij ^ represent, respectively, 
the nearest, next-nearest and next-next-nearest neighbor 
pairs, and Ji 3> J2, J3 > 0. The summation J2fj.u e{a} 
taken over sites in the a-th hexagon. The operators P4 
and Pq represent, respectively, the 4-body and 6-body 
ring exchange interactions^^, 

P4 = ^ + J2s^-S, + AG,jki, (6) 

Gijki = {Si ■ Sj){Sk ■ Si) + {Si ■ Si){Sj ■ Sk) 

- {S,-Sk){SrSi) (7) 



on the honeycomb lattice are macroscopically degenerate 
ground states of Hk- An example of the dimer covering 
states is shown in FIG. 1(c). The spin-spin correlation 
functions for these VB states exhibit long-distance expo- 
nential decay, indicating the absence of magnetic long- 
range order. Also, the result of the single-mode approx- 
imation supports the existence of a spin excitation gap 
above the spin-singlet ground states. However, because 
each dimer state breaks the spatial symmetry of the sys- 
tem, a quantum spin liquid state is not the unique ground 
state of the Klein model. 




FIG. 1: (a) Honeycomb lattice, (b) Three kinds of four- 
spin exchange processes on a hexagon, (c) An example of 
a dimer covering state. The arrows indicate the phase con- 
vention of the singlet states. An arrow from i to j represents 
Oij\0). Hexagons containing circles are covered with unflip- 
pable dimers. (d) Staggered VBC. 
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where the sum in eq. is taken over all combinations of 
(i, j, fc, Z, m, n) in three different pairs depicted in FIG. 2. 
(— 1)^ is 1 (-1) when a combination is generated by an 
even (odd) number of transpositions between different 
pairs. The summations g ^ are, respectively, taken 
over the four-spin configurations of types A, B, and C of 
the a-th hexagon, as depicted in Fig. 1(b). The operator 
V^^ projects the six spins surrounding the a-th hexagon 
(depicted as i, j, k, I, m, n in Fig. 1(a)) onto the sub- 
space with total spin S = 3. The summation ^^--j 
is taken over these six sites in the a-th hexagon. is the 
Hamiltonian of the s — 1/2 Klein model on a honeycomb 
lattice, which was studied in detail by Chayes et ali^SiSi 
The ground state space of TJk is spanned by valence bond 
(VB) states, which are formed from spin-singlet pairs of 
nearest-neighbor sites. Thus, all dimer covering states 



In the following, we show that _ffR introduces a reso- 
nance among these dimer covering states and, as a re- 
sult, selects an almost unique ground state that can be 

(9) 

expressed as a superposition of VB states, preserving the 
spatial symmetry. To this end, it is convenient to uti- 
lize the Schwinger boson representation for spin opera- 
tors^S; S+ = u\d,, S- = d\u,, = {u\u, - 4d,)/2, 
and Si = {ului + d\di)/2. Here the average value of 
S is 1/2. We also introduce the spin-singlet operator 
Oij = {uidj — diUj)/\/2. Then, each dimer state is ex- 
pressed as 



\D)= n 



(10) 



Here, T) is the set of nearest-neighbor pairs correspond- 
ing to a particular realization of the dimer covering. 
We choose the phase convention of Oij as depicted in 
FIG. 1(b). 

We now show that is a variant of the QDM repre- 
sented by spin operators. The Hamiltonian of the QDM 
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FIG. 2; The spin-singlet states of a hexagon. 



-.18.22 



ffQDM = -i[|l)(2| + |2)(1|] + y[|l>(l| + |2)(2|]. (11) 

Here, |1) and |2) are the dimer covering states of a single 
hexagon corresponding, respectively, to the states (1) and 
(2) of Fig. 2. The first term in TJqdm is the kinetic term 
that transfers |1) to |2) and |2) to |1). The second term 
is the potential term. To construct operators that act as 
the first and second terms of Hqdm and are expressed 
in terms of spin operators, we consider all possible spin 
singlet states of a single hexagon, depicted in Fig. 2. Al- 
though the states except (1) and (2) in Fig. 3 are not the 
ground state of Hk, the kinetic term for the states |1) 
and 1 2) are expressed by the linear combination of these 
states. To see this, we introduce the density operators for 
the singlet dimers on a hexagon corresponding to these 
fifteen states, 

rji) = ol^o.,ol^o,^ol^o,i, (12) 



ri^) = ol^col^Ox^ol^o^^, (13) 

and so forth. It is obvious that T('")|m) = \m) with 
m=l,2,...,15. In particular, we are concerned with the 
action of Ta^^ and Ta^^ on the states |1) and |2), i.e. 



rW|i) = |i), rj2)|2) = |2), 
1, 



1 



r(2)|i\ 



rl2>. (14) 



The first pair of relations here implies that Ta^'^ and T^^^-* 
act as the potential terms of Hqbm- In the following, it 
is convenient to use the representations 



(15) 



with 



hia = O.^OakO,,?, (16) 

h2a = COacOw (17) 
Next, we introduce the operator, 

/C, =-/lI>2a-4a/lla, (18) 

which is (as can be seen by using the relations 0^^^0\^ + 
= OlMx^ etc.) expressed as 
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Note that JCa satisfies the following relations: 

^a|l)--|2> + i|l), /Ca|2)=-|l)-Ki|2). (20) 

These actions of ICa are similar to those of the kinetic 
term of Hqum ■ To establish the complete set of relations 

between the QDM and the operators ICa, Ta^^ and Ta'^\ 
we need to verify their actions on unflippable dimers that 
are not flipped by ICa, as shown in FIG. 1(b). It is easily 



shown that the actions of IC and T^^^ + T^^ 



the 



flippable dimers generate dimer covering states that are 
not in the ground-state space of -ffx- These non-ground 
states contain at least one singlet pair in the set of six 
spins surrounding the hexagon to which an operation is 
applied; i.e. the spins on i, j, k, I, m, n in FIG. 1(a) when 
the operators are applied to the a-th hexagon. Thus, 
the unwanted states are excluded by the projection onto 
the maximum spin states of these six spins. This is 
carried out by applying the operator P*^. Also, using 
the relation Si ■ Sj = j — 0\jOij, we find ICa = Ka, 
Ta^^ + = Va, thereby arriving at the Hamiltonian 



B. Exact ground state 

We now show that the Hamiltonian H — Hk + 
has the same ground state properties as the QDM in 
a certain parameter region. The lowest energy sector in 
the Hilbert space of H is spanned by VB states, provided 
that the ground state of can be expressed as a linear 
combination oi \ D). To find the ground state, it is useful 
to rewrite Hn into the form 

Hn = '^^Y.'P'Jihl-hDihia-h.a) 



2 

J3-J2 



J2'^a'ihl+hl){h,a+h2a). (21) 
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For J3 > J2, the energy eigenvalues of H^. are non- 
negative, and only unflippable dimer states are zero en- 
ergy states. Therefore the ground state of -ff is a dimer 
covering without flippable dimers. If we simply impose 
periodic boundary conditions in both x and y direc- 
tions, the ground state remains macroscopically degen- 
erate. This degeneracy is removed by imposing periodic 
boundary conditions with a shift by one hexagon in the y 
direction. Then, the unique ground state is the staggered 
VBC (see FIG. 1(d)), which was previously identified by 
Moessner et al^ 

The energy levels of Hn arc non-negative also in the 
case of J2 = J3. Although the staggered VBC is obvi- 
ously a zero energy state here too, in this case there also 
exists a non-trivial, liquid-like ground state. It is eas- 
ily seen that the equal-amplitude superposition of dimer 
states \G) = ^ - \Di) is the zero energy state of if 
the summation is restricted within a sector of dimer 
states related by the local operations ifT^ and iPHjl . Gen- 
erally, a class of dimer states connected by these local 
operations is characterized by a topological numberiiSiSS. 
On a honeycomb lattice, however, in addition to a topo- 
logical number, the total number of dimers in each direc- 
tion is also conserved by the local flip. Thus, \G) is the 
unique ground state of in a given sector specified by 
these conserved quantities. This state, preserving both 
the spin-rotational and spatial symmetries, is a spin liq- 
uid state, and it is in the same universality class as the 
short-range RVB state. In fact, the upper limit of the 
spin-spin correlation function for this state exhibits long 
distance exponential decay expressed bji2& 



|(G|5f^||G)|<(3/2)2-l'-- 



,1 



(22) 



These results show that the QDM is realized as a spin 
system with multiple ring exchange interactions, and the 
case J2 = J3 corresponds to the Rokhsar-Kivclson (RK) 
point of the QDMi^. 

In the case of J2 > J3 , unfortunately, we have not been 
able to derive an exact ground state of H analytically. 
However, the plaquette VBC state obtained by Moessner 
et al.-- for the QDM in this parameter region is indeed 
an exact eigenstate of H, and it is possible that this is 
the ground state. 

We now consider the Hamiltonian H without the pro- 
jection operator V^^ in H^. Then, the model is more 
realistic, though its exact ground state is no longer ac- 
cessible by analytical method. The operation of Ka and 
Va on unflippable dimers do not vanish, but create states 
which are not in the ground state sector of H^. In this 
situation, the staggered VB solid state is not the zero- 
energy state even for J3 > J2, and pushed up to a state 
with higher energy of order Ji. Thus, it is expected that, 
at the RK point J2 — J3, the spin liquid state is the 
most plausible candidate for the true ground state pro- 
vided that Ji ^ J2, J3. 



C. Dimer-dimer correlation in the spin liquid state 

In contrast to the magnetic correlation, the dimer cor- 
relation in the spin liquid state obeys a power law, as 
shown below. In the ground state space, the dimer- 
dimer correlation is identical to that of the classical dimer 
model on a honeycomb lattice, which belongs to the uni- 
versality class of the Gaussian model with central charge 
c = 1 The correlation function for the dimer density 
operator at the bond (ij), Nij = —Si ■ Sj + 1/4, displays 
the long-distance behavior 



[cos(47r|r, -r;|/3)-l] 



(23) 



when the two dimers on (ij) and (Im) are in the same 
direction^*^ . This power law decay implies that the low- 
energy properties are governed by gapless non-magnetic 
excitations. The low-lying excitation energy is computed 
using the single-mode approximation. We assume that 
the excited state takes the form \k) = e^''^^6Nij\G), 
with SNij = Nij 



(Nij). By definition, |fc) is orthogo- 
nal to the ground state; i.e. {k\G) = 0. Therefore, the 
excitation energy is obtained from 



{G\[6N_k[H,6Nk]]\G) 
{G\5N.kSNk\G) ■ 



(24) 



The denominator of the right-hand side of this relation 
can be computed from the correlation function (NijNki) 
and behaves as ~ ln(i?/a), where a and R are the lat- 
tice constant and the system size, respectively. When 
the excitation energy is sufficiently smaller than the spin 
excitation gap of H^, the total number of spin-singlet 
dimers is conserved, and also, Hk does not contribute 
to the low-lying excitations. Thus, the numerator of the 
right-hand side of the above relation for Sk reduces to 
J2ii(G\SNijHnSNimr,ri\G)P for small /c. The coefRcient 
of the fc^-term is expressed in terms of the three-body 
dimer correlation functions, which exhibits a leading log- 
arithmic behavior of the form ~ ln(i?/a). Hence, the log- 
arithmic divergences in the denominator and numerator 
of efc cancel out. We thus find that the gapless excitation 
energy in the spin-singlet sector is given by Sk ~ ak^. 
This dispersion relation is in accordance with Henley's 
conjecture based upon height representations^. 
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FIG. 3: (a) Dimer covering on a square lattice, (b) A ground 
state of Hk that is not a ground state of Hr. (c) The stag- 
gered VB crystal. 
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III. s = 1/2 HEISENBERG 
ANTIFERROMAGNET ON A SQUARE LATTICE 

The analysis presented in the previous sections is easily 
extended to other lattice structures. For the model with 
s=l/2 on the square lattice, the Hamiltonian is also given 
by H = Hx + -ffR, but in this case we have the following: 



Hk = 



1 E iS.-S,){S,-S,) + ^N], (25) 



j^k,l 
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with Qa = Z]ijG{a} Si-Sj. Here the summation I]ye{a} 
is taken over the nearest neighbor sites of the a-th square 
(the eight sites on i,j ,k,l,m,n,o,p in FIG. 3(a)). The op- 
erator P°4 is the 4-body ring exchange interaction on the 
a-th square. The ground state space of the s = 1/2 Klein 
model on a square lattice is more complicated than that 
of the model on a honeycomb lattice. In the case of a 
square lattice, in addition to dimer covering states, some 
states with spin-singlet pairs on next-nearest neighbor 
sites are also in the ground state space. An example of 
such states is depicted in FIG. 3(b). However, fortunately, 
such singlet states are not in the ground state space of 
i?R, and they can be excluded from our consideration. It 
is thus found that at the RK point J2 — J3, under peri- 
odic boundary conditions in both x and y directions, the 
ground state of Hk + is the short-range RVB state 
J2d I^)' where \D) is the dimer covering state on the 
square lattice. The low energy properties of this short- 
range RVB state were investigated in detail by Rokhsar 
and Kivelsoni^. For J3 > J2, the staggered VBC state 
depicted in FIG. 3(c) is the unique ground state. 



IV. s = 1 HEISENBERG ANTIFERROMAGNET 
ON A TRIANGULAR LATTICE 

In this case, we divide the triangular lattice into 
three honeycomb lattices, hj^, and he, as shown in 
FIG. 4(a). Here, each site consists of two sites of two 
different honeycomb lattices. The Hamiltonian on this 



lattice is also given by = Hk + i?R. The Klein Hamil- 
tonian Hk for this system has the same form as that 
given in Eq.Q, but, here it is defined on a triangular 
lattice, and it is expressed in terms of the s — 1 spin 
operators. The ground state space of this Klein model is 
spanned by states fully packed with loops composed of 
singlet dimers on the three kinds of hexagons arranged 
sequentially in order the A, B, C, A, B, C,... We show 
an example of loop covering states in FIG. 4. Resonance 
among these loop covering states is introduced as. 



Hk = Ha + He 



He 



(30) 



where Ha, Hb and He are defined on Ha, Hb and he, 
respectively, and take the form of Eq.(|2Jl with re- 
placed by the projection onto the maximum spin state of 
the thirteen nearest neighbor spins of the a-th hexagon 
on the parent triangular lattice: 



V. 



M 



n 

s=0 



(E 



ie{a} 



182 - s{s + 1) 



(31) 



Here the summation Ejg{a} taken over the nearest 
neighbor sites of the a-th hexagon. (See FIG. 4(b).) 

The ground state properties are similar to the mod- 
els considered in the previous sections. At the RK point, 
the equal-amplitude superposition of loop covering states 
is the unique ground state in a given sector. The loop 
statistics of this state are described by a classical loop 
model referred to as the "red-green-blue model". The 
loop correlation functions of this state exhibit power-law 
long-distance behavior—. Therefore the excitation en- 
ergy in the spin singlet-singlet channel of this system is 
gapless. This ground state represents a new universal- 
ity class of a spin liquid described by the quantum loop 
model. 



V. DISCUSSION AND SUMMARY 

In this paper, we have presented rigorous results for 
two-dimensional quantum spin systems with multiple 
ring exchange interactions that possess quantum spin liq- 
uid ground states. It should be noted that Misguich et 
al. and LiMing et al. showed that for the s = 1/2 Heisen- 
berg antiferromagnet on a triangular lattice, sufhciently 
strong ring exchange interactions destroy magnetic long- 
range order, and stabilize a spin liquid state ^Siii. The 
spin liquid state found by them has a remarkable simi- 
larity with the ground state discussed in the present pa- 
per: exponentially decaying spin-spin correlations, and a 
large number of spin singlet excitations in the spin gap 
(at least for some parameter regions) . Since the classical 
Heisenberg model with ring exchange interactions does 
not possess macroscopic degeneracy, we speculate that 
the underlying structure of the model may have a close 
connection with the Klein model, which is the origin of 
singlet gapless excitations in our models. 
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FIG. 4: (a) Triangular lattice divided into three honeycomb 
lattices, /ia (black solid), hs (black dotted), and he (grey). 
The bold lines represent an example of loop covering, (b) 
Open circles represent the nearest neighbor sites of the a-th 
hexagon. The a-th hexagon is denoted by thick lines. 



These observations imply that the realization of the 
spin liquid state due to ring exchange interactions may be 
generic in quantum antiferromagnets, though the models 
considered in the present paper is rather complicated. 
The multiple ring exchange interaction is a result of 
strong quantum fluctuations, which destabilize conven- 
tional magnetic orders. In real systems, effects of such 
strong quantum fluctuations have been extensively stud- 
ied so far in connection with magnetism of solid '^He2^. 
Also, it has been discussed that multiple ring exchange 
interactions may affect significantly the magnetic struc- 
tures of NiS2. Thus, there is a possibility that the quan- 
tum spin liquid state due to ring exchange interactions 
may be realized in real systems with sufficiently strong 
multiple ring exchange interactions. 
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